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Abstract 



The HERMES hydrogen and deuterium nuclear-polarized gas targets have been 
in use since 1996 with the polarized electron beam of HERA at DESY to study 
the spin structure of the nucleon. Polarized atoms from a Stern-Gerlach Atomic 
Beam Source are injected into a storage cell internal to the HERA electron ring. 
Atoms diffusing from the center of the storage cell into a side tube are analyzed 
to determine the atomic fraction and the atomic polarizations. The atoms have a 
nuclear polarization, the axis of which is defined by an external magnetic holding 
field. The holding field was longitudinal during 1996-2000, and was changed to 
transverse in 2001. The design of the target is described, the method for analyzing 
the target polarization is outlined, and the performance of the target in the various 
running periods is presented. 



1 Introduction 



The HERMES experiment at HERA (DESY, Hamburg) was designed to study 
the quark-gluon spin structure of the nucleon by deeply inelastic scattering 
(DIS) of longitudinally polarized electrons from an polarized gas target of 
hydrogen and deuterium atoms internal to a storage ring[l,2]. The 27.5 GeV 
electron beam stored in the HERA storage ring with currents up to 50 mA 
is polarized by the Sokolov-Ternov effect [3], with a time constant of about 
30 minutes. A longitudinal spin direction is obtained where the experiment is 
located within the HERA east straight section by means of a pair of spin rota- 
tors [4]. In a machine with one pair of spin rotators before the 2001 upgrade, 
a typical polarization of 60% was achieved. 

The HERMES experiment [5] was proposed in 1990, and is designed to run 
in parallel with the ZEUS and HI collider experiments without significantly 
reducing the electron beam lifetime, which is of the order of 10 hours. This 
limits the target areal density to about 10^^ H-atoms/cm^ and excludes the 
use of any solid target material. A target employing a cold storage cell fed by a 
a polarized atomic beam source was proposed and implemented. Polarized gas 
targets for storage rings were reviewed recently by Steffens and Haeberli [6]. 

Since 1996, the present target apparatus for polarized hydrogen and deuterium 
has been in operation. During the years 1996-97 longitudinally polarized hy- 
drogen was used, while longitudinally polarized deuterium was employed 
in 1998-2000 and from 2002 onwards, transversely polarized hydrogen 

{H±) has been in use. 

^ Corresponding author: Tel:-|-49-40-89981963; e-mail: Paolo.Lenisa@desy.de 
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Apart from a high areal density and rehable operation over many years, the 
most challenging requirement for the internal target was to determine the ef- 
fective target polarization seen by the beam. The target polarization cannot 
be measured by making use of scattering processes involving the high-energy 
electron beam due to the small cross sections involved. Nuclear vector po- 
larization Pz and, for deuterium also the tensor polarization P^^ were to be 
measured, as well as the electron polarization Pe for diagnostic purposes. The 
required accuracy was 5P = 0.03 absolute uncertainty on the average po- 
larization over long measuring periods, and 1 % statistical accuracy within 
a couple of minutes in order to monitor possible changes in time. This goal 
was achieved by employing a sampling technique in conjunction with detailed 
simulations and frequent calibration measurements. 

In the present article, an overview of the HERMES polarized hydrogen and 
deuterium gas target is given with reference to the various papers describ- 
ing the subsystems of the target. The atomic beam source (ABS) [8], the 
storage cell [9], the Target Gas Analyzer (TGA) [10] and the Breit-Rabi po- 
larimeter (BRP) [11] are summarized in Section 2. Section 3 describes the 
analysis applied to the raw TGA and BRP signals and the several calibra- 
tions required to derive measured values for the degree of dissociation and 
the atomic polarizations of the sample gas ^ . Section 4 reviews the various 
physical processes inside the storage cell that directly influence the TGA and 
BRP analyses [12,13,14,15]. The results of the target analysis are reported in 
Sections 6 and 7. Section 8 mentions a selection of measurements carried out 
by HERMES that crucially depended on the polarized target. The paper is 
concluded in Section 9. 



2 Survey of experimental apparatus 

The HERMES polarized hydrogen and deuterium target, schematically shown 
in Fig. 1, consists of an atomic beam source (ABS) that injects a spin polarized 
beam of H or D atoms into a thin-walled storage cell. A sample of gas diffuses 
from the center of the cell into a Breit Rabi Polarimeter (BRP) that measures 
the atomic polarizations, or into a Target Gas Analyzer (TGA) that measures 
the relative atomic and molecular content of the gas. A magnet surrounding 
the storage cell provides a holding field defining the polarization axis and 
preventing spin relaxation by effectively decoupling the magnetic moments of 
electrons and nucleons. 



^ In this article, the term "atomic polarizations" refers to the polarizations of the 
electrons and nuclei in the atoms 
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Fig. 1. Schematic view of the HERMES longitudinally polarized target. From left 
to right: Atomic Beam Source (ABS), target chamber with cell and magnet, and 
diagnostic system composed by Target Gas Analyzer (TGA) and Breit-Rabi Po- 
larimeter (BRP). The locations of the radio- frequency transition (RFT) units are 
indicated. 

2.1 Target chamber 



The target chamber shown in Fig. 2 is evacuated by two turbo-molecular 
pumps with a combined speed of 4400 ls~^. During operation the chamber 
pressure is typically in the 10~^mbar range due to the high gas load from 
the atomic beam, which enters the chamber from 30° above the horizontal 
plane. The storage cell and its supporting cooling rails are cantilevered from 
the upstream^ end flange. 

A 0.3 mm thick stainless steel exit window on the downstream end of the tar- 
get chamber allows the scattered electrons and hadronic fragments to leave 
the target chamber and to be detected by the HERMES spectrometer [5]. This 
thin window is integral with the following thin-walled beam tube with no in- 
tervening flange, in order to minimize the high-Z material in this region and 
reduce shower production by 7r° mesons produced at small angles. About 1 m 
upstream of the target cell, an adjustable collimator (labeled CI) is located to 
stop direct synchrotron radiation. It normally has the smallest scaled aperture 
in the HERA ring, and is opened during electron injection [5]. A fixed collima- 
tor (labeled C2) just upstream of the storage cell protects it from synchrotron 
light scattered by CI and from exposure to electron showers from CI. 



Upstream and downstream mean with respect to the HERA lepton beam direction 
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Fig. 2. Longitudinal target chamber and superconducting magnet viewed from 
downstream with respect to the HERA beam direction (left), and from above (right). 

For longitudinal running the storage cell and vacuum chamber are immersed 
in a longitudinal field generated by a superconducting magnet. The magnet 
consists of four coils and an iron yoke designed to provide optimum field uni- 
formity over the length of the storage cell while providing a central gap for the 
atomic beam to enter the cell. The magnetic field (up to 350 mT) provides the 
quantization axis for the spins of the polarized atoms in the storage cell, and 
decouples the spins of the nucleons and electrons. The passing HERA beam is 
bunched to provide high luminosity for the collider experiments HI and ZEUS. 
As discussed in Sec. 4.2.4, this may induce resonant nuclear depolarization of 
the target gas. The induced resonances can be avoided by proper setting of 
the holding field, provided that its homogeneity is sufficient [13]. In the case 
of hydrogen, the maximum permissible variation of the magnetic field along 
the storage cell is 1.5%. For operating with deuterium gas, this requirement 
was found to be less critical [16]. 



2.2 Transverse magnet 



After the end of the longitudinal running in 2000, the HERMES target was 
modified for measurements with a hydrogen gas target in a transversely ori- 
ented magnetic holding field. The intensity of the transverse magnetic field 
was limited by the amount of synchrotron radiation power generated by the 
defiection of the beam by the target magnet (5kW maximum), so it was de- 
cided to build a normal conducting magnet. As will be explained in Sec. 4.2.4, 
a better homogeneity of the transverse target magnet of AB < 0.15 mT at a 
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Fig. 3. Upper panel: field strength of tlie longitudinal target magnet measured along 
the HERA beam direction z, covering the full length of the storage cell. The line 
shows a polynomial fit to the measurement. Lower panel: transverse magnet field 
uniformity measured along the z axis at the nominal field strength B = 297 mT. De- 
viations of ABy = 0.15 mT (vert.) and AB^ = 0.60 mT (horiz.) have been measured 
within the cell volume. 

field value of about 300 mT was required. 



The geometrical constraints imposed by the HERMES acceptance made it 
difficult to design a magnet fulfilling the uniformity requirement. After con- 
struction, the field uniformity was measured and improved by shimming the 
pole tips. At a field strength of 5=297 mT, maximum deviations of AB^ = 
0.05 mT, ABy = 0.15 mT and AB^ = 0.60 mT within the storage cell volume 
have been achieved, z being the longitudinal direction, y the vertical and x 
the horizontal axis. The magnet was installed in July 2001. The strength of 
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period of use 


cross section 

Ui /\ U I llilli I 


conductance 


05.96 - 11.99 


29.8 X 9.8 


lA5y^T/M 


01.99 (test) 


19.8 X 8.0 


0.61^yT/AI 


12.99 - now 


21.0 X 8.9 





Table 1 

The history of cells used by HERMES. The lengths of the major and minor axes of 
the elUpse are indicated by a and 6; T is the temperature of the cell in K, and M 
is the atomic mass of the gas in amu. 

the transverse field vs. z is shown in the lower plot of Fig. 3. 

2.3 The storage cell 

The HERMES storage cell [9] is made of two 75 /xm thin pure aluminum sheets, 
which are tightly spot- welded together. It is 400 mm long and has an elliptical 
cross-sectioned shape of 21 x 8.9 mm^. All the aluminum parts of the cell are 
coated with Drifilm [17] to minimize gas depolarization and recombination 
caused by wall collisions [18,19,20,21]. A vented extension of the cell contin- 
ues about 147 mm downstream of the cell before joining the thicker- walled 
wakefield-suppressing transition tube, to ensure that all scattered particles 
inside the spectrometer detector acceptance pass through the minimum thick- 
ness of material (see Fig. 4). The wake field suppressors are made of 100 /xm 
thick titanium, and confine the electromagnetic high frequency field of the 
bunched HERA beam, forming a gradual electrical transition between the 
storage cell and the HERA beam pipe. 

Cells with three different elliptical cross sections have been used thus far. The 
mentioned size of 21 x 8.9 mm^ is the one in use since December 1999, called 
the medium-sized cell. The largest was used until December 1999, and had 
a cross section of 29 x 9.8 mm^. (A small one (19 x 8.0 mm^) was tested in 
January 1999 and found to be too small.) Tab. 1 summarizes dimensions and 
conductances of the three cells that have been used. The conductance values 
apply to gas injected into the cell center. 

The storage cell is mounted onto cooling rails, which are cooled by temperature 
regulated gaseous helium. The optimum operating temperature for hydrogen 
targets was found to be about 100 K, where recombination and depolarization 
effects are low. Additionally, the lower conductance at this temperature leads 
to a higher target thickness by a factor -\/3 compared to room temperature. 
For deuterium gas the target temperature could be decreased further to about 
65 K, which increased the target thickness by about 23 % compared to 100 K. 
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Fig. 4. The storage cell and its support flange. 



The temperature of the cooUng rails is monitored at three locations by plat- 
inum resistors. Since 1999 an additional thermocouple wire has been attached 
to the aluminum foil close to the cell ellipse to directly monitor the cell tem- 
perature during injection of the HERA beam when a temperature rise of about 
5-10 K is usually detected. The cooling rails are mounted on the target cham- 
ber flange on the upstream end. Due to thermal contraction of the cell when 
being cooled, the position of the cell center changes relative to the atomic 
beam source and the Breit-Rabi polarimeter. From room temperature down 
to 100 K the cell center moves by 0.9 mm. The alignment of the atomic beam 
source and Breit-Rabi polarimeter relative to the cell center includes an offset 
to compensate for this effect at the target operating temperature. 



The cell is shown along with its support flange in Fig. 4. Two side tubes are 
connected to the beam tube, one to inject the polarized atoms and another 
smaller one to sample about 5% of the target gas and analyze it in the target 
gas analyzer and the Breit-Rabi polarimeter. The sample tube is tilted by 
120° with respect to the axis of the injection tube to ensure that the sampled 
atoms have thermalized with the storage cell wall. A permanent dipole magnet 
around the sampling tube outside the target chamber extends the magnetic 
holding field of the target magnet to prevent depolarization (see Fig. 2). A 
capillary near the location of the feed tube allows injection of unpolarized gas 
into the cell. 
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Fig. 5. Schematic view (from downstream of the electron beam) of the HERMES 
ABS with dissociator and colhmator for beam formation. Two sets of sextupole 
magnets are located along the beam axis as are the high-frequency transitions. The 
axis of the ABS is tilted by 30° downwards with respect to the horizontal plane. 

2.4 Atomic Beam Source 



The atomic beam source [8] consists of a dissociator, a powerful differential 
pumping system, a beam forming system, a sextupole magnet system to focus 
atoms with = +i into the storage cell, and adiabatic high-frequency tran- 
sitions to manipulate the hyperfine population of the atomic beam. Injected 
fluxes of ^^^^ ^ 6.5 X 10^^ atoms/s in the case of hydrogen (2 states) and 
^ABS ^ ^ j^gie atoms/s in the case of deuterium (3 states) have been ob- 
served. These values are calculated from the density- dependent spin exchange 
relaxation [16]. The injected nuclear polarization P^"-' was constantly above 
0.97 for hydrogen and above 0.91 for deuterium. A schematic diagram of the 
ABS is depicted in Fig. 5. 

Molecular hydrogen/deuterium gas enters the dissociator from the left. The 
molecules are dissociated by a radio frequency discharge with a frequency of 
13.56 MHz in a pyrex tube, producing a degree of dissociation up to 80 % (H2, 
D2) at a flow of about 1 mbarls"^ and a radio frequency power of 300 W. To 
maximize the stability and the degree of dissociation, an amount of oxygen 
between 0.1 and 0.3 volume percent is added to the molecular gas. 

A high degree of dissociation at higher flow can potentially be achieved using 
a recently developed microwave (MW) dissociator [22], which operates at a 
frequency of 2.45 GHz. With typical flows between 1 and 2mbarls~^ and a 
microwave power of about 600 W, a degree of dissociation higher than 80 % 
is achieved. The MW dissociator was used in the year 2000 when running 
with deuterium, when it provided a target thickness improvement of 5-10 %. 
The RF dissociator has been used since, because the ABS sextupole system is 
not optimized for the lower atomic velocity distribution produced by the MW 
dissociator, eliminating the possible benefit for hydrogen. 



11 



The atomic gas flows through a conical nozzle with an opening diameter of 
2 mm, which is cooled to about 100 K. The water produced in the discharge 
freezes in the nozzle. Although this beneflcially reduces surface recombination, 
the ice layer continuously thickens and slowly reduces the gas conductance of 
the nozzle. After typically 3-5 days of operation, the pressure in the disso- 
ciator tube rises substantially, increasing volume recombination in the tube 
and resulting in a noticeable decrease in the degree of dissociation of gas ex- 
iting the nozzle. Full recovery is achieved by slightly warming up the nozzle 
to remove the ice layer. 

After passing the cold nozzle, the gas expands into the vacuum of the disso- 
ciator chamber. A powerful differential pumping system with a total nominal 
pumping speed of more than 15000 ls~^ suppresses scattering of the atomic 
beam, and ensures low gas flow into the HERA vacuum system. 

The magnet system of the HERMES atomic beam source consists of a total of 
flve Halbach type [23] segmented permanent sextupole magnets with a maxi- 
mum poletip fleld of about 1.5 T. Every magnet consists of 24 segments of Va- 
codym alloy ^ . Four high-frequency transition units are available: two between 
the sextupole magnets and two after the last magnet. The system provides 
enough flexibility that many different combinations of hyperflne states can be 
injected into the cell. Between the sextupoles, a strong fleld transition (SET) 
and a medium fleld transition (MET) can be used simultaneously, whereas 
the weak fleld transition (WET) and the SET after the sextupole magnets 
can only be used separately, because they share the coils for the static and 
gradient magnetic flelds. The complete set of possible hyperflne populations 
that can be produced is shown in Tab. 2. 



2.5 Target Gas Analyzer 

The target gas analyzer (TGA) [10] measures the atomic and molecular con- 
tent of the gas extracted from the storage cell through the sample tube. The 
TGA arrangement, shown in Eig. 7, consists of a pair of baffles, a chopper, 
a 90° off-axis quadrupole mass spectrometer (QMS) with a cross beam ion- 
izer and a channel electron multiplier (GEM) for single ion detection. The 
TGA is integrated into the vacuum system of the BRP sextupole chamber, 
which is served by two cryopumps and a titanium sublimation pump, with 
a total pumping speed of about 7000 Is^^. During operation, the pressure in 
the TGA detector is about 4-10~^mbar. Prior to normal operation, the TGA 
vacuum chamber is baked with temperatures up to 180 °G for 48 hours. The 
BRP/TGA vacuum scheme is shown in Eig. 8. 

^ Vacodym is a registered trademark of Vacuumschmelze GmbH, Hanau, Germany 
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Table 2 

Injection modes of the atomic beam source. The table shows the high frequency 
transition (HFT) units employed (2"^^ and 3'''^ columns), the hyperfine states in- 
jected in the target cell and the resulting electron (Pe), nuclear (P^) and tensor 
{Pzz) polarizations in the ideal case of 100% efficiency of the sextupole system and 
transition units, no depolarization inside the target cell and infinitely strong guide 
field. The last column clarifies the purpose of each mode (Cal: calibration, moni- 
toring; Data: normal HERMES operation). The injection modes for hydrogen and 
deuterium single states are listed separately. For some modes the MFT between the 
sextupoles has to be tuned as a WFT 1-4. The no state injection mode is possible 
only for hydrogen and has been operational since April 2002. The adopted labeling 
for the hyperfine states is explained in the Breit-Rabi diagram of Fig. 6. 

The TGA is mounted 7° off-axis with respect to the BRP, in order to not 
interfere with the beam entering the polarimeter. The first baffle ensures that 
only gas from the sample beam can reach the detector, while the second one 
acts as a collimator, so that no atoms can hit the metal parts of the ionizer 
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Fig. 6. Hyperfine energy levels of hydrogen (left) and deuterium (right) atoms as 
a function of the magnetic holding field (Breit-Rabi diagram) and corresponding 
labeling. The field values are scaled with the corresponding critical field and energy 
values with the corresponding hyperfine energy. 



TGA system ^ CEM 




BRP 



Fig. 7. A schematic diagram of the Target Gas Analyzer. The gas of atoms and 
molecules diffuses through the extension tube and is collimated by two baffles before 
entering the ionizing volume. The chopper is used for background subtraction 

and recombine. The chopper, placed between the two baffles, rotates at a 
frequency of 5.5 Hz, thus periodically blocking the sample beam in order to 
allow subtraction of the residual gas signal. Particles entering the detector are 
ionized by 70 eV electrons, mass filtered with the QMS, and finally detected 
by the CEM. The pulses are counted by a 2048-bin time-resolving counter 
(TRC). The bin length is typically set to 175 fis so that two chopper turns are 
collected within the whole bin range as shown in Fig. 9. The synchronization 
between chopper status and counting rate is given by a trigger signal coming 
from the chopper motor. The beam rate is calculated by taking the difference 
between the number of counts collected with the chopper open and closed. 
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Fig. 8. Schematic view of the BRP/TGA vacuum system. 
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Fig. 9. Time spectrum of atomic count rates in the TGA. The vertical hues indicate 
the binning used to define the chopper positions open, undefined and closed. 

2.6 Breit-Rabi Polarimeter 



The Breit-Rabi polarimeter (BRP) [11] measures the relative populations Ui 
of the hyperfine states of hydrogen (or deuterium) atoms contained in the 
sample beam. From this measurement, the absolute atomic polarizations can 
be calculated by applying the knowledge of the target magnetic field strength. 

A schematic view of the BRP is shown in Fig. 10. From left to right, the 
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Fig. 10. Schematic layout of the BRP. The rf-transition units extend from z = 600 
to 810 mm. The beam blocker at the entrance to the first 6-pole magnet ensures 
100% rejection of atoms with = — ^- The beam shutter is used to measure the 
hydrogen contribution coming from dissociative water ionization. 

sample beam leaves the extension tube of the target cell, encountering first 
two hyperfine transition units, then a sextuple magnet system and eventually 
the detector stage. During operation, a differential pumping system keeps the 
pressure at 2xl0~^mbar in the sextupole chamber and at 2.5xl0~^° mbar in 
the detector chamber. 



Two transition units are used to exchange the populations between pairs of 
hyperfine states: a strong field transition unit (SFT) with tilted resonator 
that can be tuned for both tt and a transitions, and a medium field transition 
(MFT) unit that can induce various vr transitions according to the static field 
strength and gradient setting used. The sextupole system is composed of two 
magnets which focus the atoms of the sample beam with = +| towards the 
BRP geometrical axis and defocus the ones with = — |. A 9 mm diameter 
beam blocker placed in front of the first sextupole magnet ensures that no 
atoms in = — | states can reach the detector. The sextupole system was 
redesigned in the year 2000 in order to improve the transmission probability 
at the typical cell temperature. The detector stage is identical to the one 
employed for the target gas analyzer: a cross beam ionizer, a quadrupole mass 
spectrometer (QMS) and a channel electron multiplier (CEM). In contrast 
to the TGA, only hydrogen (or deuterium) atoms are detected by the BRP. 
The chamber residual gas correction is carried out by using a chopper in the 
same way as described for the TGA. An additional beam shutter is placed just 
after the extension tube exit to measure the possible contribution of atomic 
hydrogen coming from the dissociation of H2O and H2 (D2O, DHO, D2 and 
DH in the deuterium case) present on the vacuum chamber walls ^ . 



^ The beam shutter was installed at the beginning of the deuterium running period 
in 1998. During the 1996/'97 hydrogen running period, the contribution of dissocia- 
tive hydrogen ionization was taken into account by detecting the intensity of masses 
2 and IBamu in addition to mass 1 amu [16]. 
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3 Analysis of TGA and BRP signals 



3. 1 Interpretation of the TGA signals 



The measured TGA atomic (molecular) signal Sa (Sm) is proportional to the 
product of the particle flux $a ($m) into the TGA, the mean inverse velocity 
of the particles {^^^ ((ir~))' ionization cross section 0"^°" (o"m") ^^d the 
detection probability ea (cm)- A dependence of the signal on the emission 
current is measured and taken into account by the normalization function 
f{Iemiss)- Taking all parameters into account one finds: 

\ ^a,m I 

The total fiux ^^"^"^ of H or D nucleons into the TGA can be related to the 
particle fiuxes $a and $m by: 

$TGA ^ ^ 2$^. (2) 

The TGA geometrical acceptance is such that detected particles have bounced 
typically more than 100 times on the extension tube inner surface [16]. For 
this reason, particles entering the TGA detector have thermalized with the 
extension tube wall, which is a few tens of degrees K warmer than the cell. 
For thermalized atoms and molecules one has: 

\Vml \Va/ 

Using eqs. (1), (2) and (3), it can be shown that the partial beam rates Sa and 
Sm can be combined to form the variable Stot which is proportional to <l>^'^^: 

Stot = Sa + kV2S^ oc (4) 

where k is the calibration constant of the TGA giving the relative sensitivity 
of the TGA for atomic and molecular beams: 

n = ^^^^. (5) 



Eq. (4) applies to the case of constant extension tube temperature T^^t and 
emission current lemiss- Taking into account the dependence of Sa and Sm on 
the measured^ Text and lemiss, one can define a temperature and emission 
normalized nucleon fiow rate (ptot into the TGA ionizer, which is proportional 

^ A thermocouple on the extension tube was first installed at the beginning of 
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to the injected ABS intensity (neglecting the small contribution from 

the target chamber rest gas): 



Kt = oc . (6) 

/ \ emiss ) 

The quantity (ptot is a good monitor of ABS stabihty over time, as it is in- 
dependent of recombination inside the storage cell. Following eq. (6) one can 
also express normalized nucleon flow rates for atoms, (pa, and for molecules, 
0m, by: 



f(I ■ ) 

J y^emiss) 

\/2KSmVTe 



(7) 



ext 



Y-m f(I . ) ■ 

J y-'-emissJ 

The degree of dissociation of the sample beam, defined as the fraction of 
nucleons in atoms relative to all nucleons entering the TGA, is given by 



0a + 0m Sa + V^K Sm 

The degree of dissociation a-^"-^^, also called atomic fraction, is displayed by 
the TGA roughly once per minute. 



3.2 The calibration of the TGA 



The purpose of the TGA calibration is to determine the relative sensitivity k, 
defined by eq. (5). The calibration can be performed by varying the relative 
amount of atoms and molecules at a constant injected ABS total flux: 

4>tot = 0a + 0m = const. (10) 

Using the relations (7) and (8), one can write 

qI SmVText \ 

X^fi^emiss) J 1 (11) 

where the minus sign reflects the fact that an increase in 0^ must correspond 
to a decrease of 0^. Hence a fit of the normalized molecular signal 0^ as 

1999. In the 1997 analysis, the extension tube temperature was related to the cell 
temperature using the expression Text{T) = Tq + T, where Tq is a constant 

estimated from calibration studies, and T is the measured cell temperature in K. 
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target /year 


K 


//||(1997) 


0.551 ±0.026 


L>l|(2000) 


0.861 ± 0.024 


i/_L(2002) 


0.55 ±0.05 



Table 3 

Summary of the values of TGA calibration constant k for the years 1997-2002. 
function of the normalized atomic one (pa results in a straight line with a slope 



Because the recombination process in the HERMES target is strongly temper- 
ature dependent, a controlled variation of the atomic fraction can be obtained 
by changing the temperature of the storage cell. A second possibility is to 
take advantage of an accidental beam loss in the region of the HERMES tar- 
get which might temporarily change the surface property of the cell and cause 
a sharp drop of a'^'-^^ followed by a slow recovery. The measured values for k 
are reported in Tab. 3. 

A detailed description of the TGA calibration and its performance can be 
found in Ref. [10]. 



3.3 Interpretation of the BRP signals 



The atomic signal detected by the BRP depends on the hyperfine population 
in the sample beam, the transmission probabilities of the sextupole system and 
on the setting of the transition units. For a stable storage cell temperature, 
the transmission of the magnetic system is fixed. The transition units can be 
switched on and off, changing the intensity of the beam reaching the detector. 
The BRP signal can be described by the following expression: 

S^ = cY.{i:^^T:^h. (12) 

where is the sample beam intensity (in atoms/s) of the atoms in hyperfine 
state I a), T^j, describes the exchange of hyperfine states by the transition 
unit z, ah is the sextupole transmission probability for state | 6), and c is a 
proportionality factor ® . The goal of these measurements is to determine the 
four (six) sample beam intensities, Jq, for the hydrogen (deuterium) sample 

If two units are working, the matrix T* has to be replaced by T^;^ = ^ac'^cfo 
^ The expression (12) simplifies into Soff = (^ala in the case that no transition 
unit is switched on. 
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beam. Summing eq. (12) over index b one obtains: 



S, = cY,M,J,, (13) 

a 

where the measurement matrix Mia includes information concerning transmis- 
sion probabihties and RF-transition efficiencies. 

Assuming that Mia is known, one can solve the system of equations (13) for la 
if, during a given ABS injection status, the BRP transition units are operated 
in at least four (six) different modes in the case of hydrogen (deuterium). The 
solution is given by 

Ia = C~'J2Ra^Si, (14) 

i 

where ^ R = M~^. The total sample beam intensity, Itot, is obtained by the 
taking the sum over all hyperfine states: 

Itot = J2la. (15) 
o 

Finally, the relative hyperfine populations Ua are given by normalizing the 
intensities la- 

Ua = (16) 

The polarizations of the atoms in the sample beam can now be calculated from 
the target magnetic field strength, as expressed by the following vectorial forms 
for hydrogen and deuterium, respectively: 



{Pe,P,f = M^ini,n2,n3,n^f 
(Pe, P^, Pzzf = Mp{ni,n2, ns, Ui, n^, n^f , 

where the matrices Mp and Mj? contain the target holding field strength, e.g. 
for hydrogen: 

/ — 1 a 1 —a \ 
M# = , (17) 

\ 1 —a —1 a J 
where a = k = B/B^, B^ = 50.7 niT. 

The statistical uncertainties of the measured signals Si are represented by 
the diagonal covariance matrix C5. Hence, the covariance matrix Cj of the 
intensities la is given by: 

-1 



Ci = RCsR^ = [m^Cs^m) . (1 



^ If the number of collected signals is larger than the number of unknowns, one has 
to use the pseudo-inverse matrix R = {M^ G s M)'^ M'^ G s , where Gs is the inverse 
covariance matrix of the measured signals. This procedure is equivalent to a fit of 
the intensities la to the signals Si. 
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The sum over all the elements of Cj gives the statistical uncertainty Altot of 
the total intensity Itot- 

{AItotf = T.(Ci)ab. (19) 

ab 

The covariance matrix Cn of the hyperfine populations Ua is 

--5(Sf)-(S:). 

while the covariance matrix of the polarization vector can be calculated in the 
following way: 

Cp = M^'^C„(M^'^)^. (21) 

From Cp one can extract the statistical uncertainty of the polarization mea- 
surement . 

During normal running conditions, the entire automated procedure to measure 
and calculate the polarization vectors and their uncertainties took roughly 60 s 
(90s) for hydrogen (deuterium), resulting in a statistical uncertainty of about 
±0.013 for hydrogen. For deuterium the statistical uncertainty was about 
±0.031 for positive and ±0.022 for negative nuclear polarization. After an 
optimization of the BRP sextupole system, and an increase in the diameter of 
the sample tube, the statistical uncertainty of a 60 s measurement for hydrogen 
was reduced to less than 0.5 %. 



3.4 The BRP calibration 



In the previous section it was shown that the atomic polarizations can be 
calculated from the measured BRP signals assuming the knowledge of the 
measurement matrix Mia- This includes the efficiencies of the transition units 
and the transmission probabilities of the sextupole system for various hyperfine 
states. The procedure used to determine these parameters is called the BRP 
calibration. 

The basic idea of the calibration measurement is to make the system of equa- 
tions (13) overdetermined so that additional parameters like efficiencies or 
sextuple transmissions can be extracted through a fit to the data. Due to the 
fact that the transition units can be combined in several possible ways, the 
number of signals which can be collected during a polarization measurement 
exceeds the number of hyperfine states Nhfg. However, the number of un- 
knowns that can be determined in this way is still less than the total number 
of parameters. The problem can be overcome by cycling this measurement 
through several (total number Nabs) injection modes (see Tab. 2) so that an 
adequate number of signals can be collected. 



21 



BRP-Signal 


Name 


MFT-RF 


MFT 


SFT-RF 


SFT 


(Jrb 


£r 

on 


Ull 


Ol 1 


Ol 1 


Ol 1 


SFT 1-4 


sl4 


OFF 


OFF 


ON 


1-4 


SFT 2-4 


s24 


OFF 


OFF 


ON 


2-4 


MFT 1-3 


ml3cl4 


ON 


1-3 


OFF 


1-4 


Ml 1 1-3 


ml3c24 


AT 

ON 


1-3 


Ol 1 


2-4 


Ml 1 2-6 


m23cl4 


ON 


2-6 


Ol 1 


1-4 


Ml 1 2-3 


m23c24 


ON 


2-3 


Ol 1 


2-4 


MFT 1-3/SFT 1-4 


ml3sl4 


ON 


1-3 


ON 


1-4 


MFT 1-3/SFT 2-4 


ml3s24 


ON 


1-3 


ON 


2-4 


MFT 2-3/SFT 1-4 


m23sl4 


ON 


2-3 


ON 


1-4 


MFT 2-3/SFT 2-4 


m23s24 


ON 


2-3 


ON 


2-4 



Table 4 



Possible switching modes of the BRP transition units for hydrogen running. During 
a calibration measurement all signals are collected, while for a polarization mea- 
surement only a subset of them is chosen. 

3.4-1 The transition efficiencies 

We first discuss the method to determine the transition efficiencies, assuming 
a priori knowledge of the sextupole transmissions (T„ for each hyperfine state. 
Calibration of those transmissions will then be discussed in Sec. 3.4.2. 

The BRP contains two transition units, a strong field (SFT) and a medium 
field (MFT) unit. The SFT can drive two transitions for hydrogen (| 1 ) ^ | 4 ) 
and I 2 ) ^ | 4 )) and five for deuterium. The efficiencies of the SFT transitions 
in the case of hydrogen are called £^14 and £^24- The MFT is a multi-step 
transition unit which can drive two sub-transitions for hydrogen (| 1 ) | 2 ) 
and I 2 ) ^ 13)^*^) and three for deuterium. For hydrogen the MFT is tuned 
either as a | 1 ) ^ | 3 ) (which is achieved by applying the 1 1 ) ^ | 2 ) and 
I 2 ) ^ I 3 ) transitions sequentially) or as a | 2 ) ^ | 3 ) only. In the latter 
case the MFT unit is tuned to avoid the 1 1 ) ^ | 2 ) sub-transition. Although 
nearly zero, the residual 1 1 ) ^ | 2 ) efficiency for the | 2 ) ^ | 3 ) tune of the 
MFT has to be taken into account. 

Due to the limited space available for the BRP equipment, there is cross-talk 
between the two BRP transition units when both static fields are switched 



" The order of the sub-transitions depends on the sign of the gradient field, which 
is negative in the HERMES BRP [16]. 
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on, so that the tuning of the MFT is influenced by the SFT field, whereas the 
SFT is not influenced by the (low) MFT field. The MFT efficiencies have to be 
measured separately for each SFT setting, as they depend on the exact static 
field. For hydrogen this leads to an increase in the number of MFT efficiencies 
from four to eight. The various MFT efficiencies are identified by means of 
the subscripts cl4 (for SFT 1-4 static field on) or c24 (SFT 2-4 static field on) 
and by the subscripts rl3 or r23 which indicate the particular sub-transition 
being considered. For example, the notation £rrai3r23c24 represents the efficiency 
of the transition MFT 1-3, step 2-3, tuned together with the static field of the 
SFT 2-4. The total number of efficiencies N^ff is 10 for hydrogen and 41 in 
the case of deuterium. The number Nbrp of possible BRP signals related to 
various transition states is 11 for hydrogen (as shown in Tab. 4) and 29 for 
deuterium. 

If the ABS is operated in Nabs different injection modes, and the BRP tran- 
sition units switch through all possible Nbrp states, the total number of col- 
lected signals is Ns = Nabs Nbrp- In this way, recalling eq. (13), one obtains 
a system of A^^ equations containing Nu = N^jf + Nabs ^hfs unknowns , 
where N^fg = 4 (6) in the case of hydrogen (deuterium). The system can be 
solved if Nf = Ng — N^ > 0, where Nf counts the degrees of freedom and is 
given by 



In such both the hyperfine state intensities la and the transition ef- 

ficiencies e can be determined. For hydrogen, the calibration measurement 
makes use of 7 ABS injection modes, resulting in Nf = 39. In the case of 
deuterium two possible options have been developed, one with 5 ABS modes 
(Nf = 74) and the other with 6 {Nf = 97). 

The solution of the system is carried out by applying a Runge-Kutta method 
as explained in [11]. The resulting covariance matrix has the form 



where Ce contains uncertainties and covariances of the efficiencies and is used 
to calculate the systematic uncertainty of the polarization measurement. The 
resulting systematic uncertainty induced by the determination of the transi- 
tion efficiencies of the order of 1.5%. Tab. 5 summarizes the transition effi- 



During the evaluation of the transition efficiencies, the sextuple system transmis- 
sions probabilities are assumed to be constant. 



Nf = Ng-N, 
= {Nbrp - 



u 



Nhfs) Nabs - Neff . 



(22) 




(23) 
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Name 


MFT 


SFT 


Efficiency (%) 


uncertainty (%) 
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no Q 

yy.o 


i.o4 


£'s24 




9-4 




1 03 


eml3rl2cl4 




i-4 


iUi.U 


1 HQ 

i.Uo 


£ml3r23cl4 








1 41 


em23rl2cl4 


1-6 


i-4 


1 n 
-l.U 


1 1 Q 

i.lo 


£m23r23cl4 






QzL Q 


1 9"^ 


eml3rl2c24 


i-o 


9 A 
Z-4 


QQ 

yy.y 


U.4Z 


eml3r23c24 






90.3 


0.82 


em23rl2c24 


2-3 


2-4 


0.0 


1.10 


em23r23c24 






97.4 


1.04 



Table 5 

BRP transition unit efficiencies measured during the 1997 hydrogen running fixing 
the sextupole transmission ratio to r2i = 1.029. Each efficiency is shown along with 
an explanation of the MFT and SFT static field setting. 

ciency measurements performed in 1997 with hydrogen. The efficiencies of the 
transitions for deuterium running can be found in Ref. [11]. 

3.4-2 The sextupole transmissions 

In the previous section it was shown how the efficiencies Ei of the BRP transi- 
tion can be measured, assuming a priori knowledge of the sextupole transmis- 
sions CT„ for each hyperfine state. The attempt to fit both Si and all of the (T„ 
within the same algorithm fails. Nevertheless, the transmission probabilities 
can be determined by separating the fitting procedure into two steps, as shown 
below. 

In the BRP (Fig. 10), atoms with ms = —\ are either defocused or stopped 
by the beam blocker, having therefore zero transmission probability (o"3 = 
(74 = for hydrogen and (T4 = (T5 = ere = for deuterium). Moreover, due 
to the normalization (eq. (16)), only ratios of transmission probabilities are 
required for measuring the hyperfine populations. Thus the remaining required 
parameters are o"2/o"i for hydrogen and o"2/o"i and a^/ai for deuterium. 

The Stern-Gerlach force of the sextupole magnets is weaker for hydrogen atoms 
in state | 2 ) compared to state 1 1 ) for magnetic fields B of the order or less 
than the hydrogen critical field B^. = 50.7 mT, near the axis of the sextupole 
system. Thus near-axis atoms in hyperfine states 1 1 ) and | 2 ) will experience 
a different transmission efficiency and the transmission ratio ri2 = o\ja2 is 
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expected to exceed unity by a small amount. Due to the smaller hyperfine 
coupling of deuterium, B^. = 11.7 mT, the deviation from unity of the ratio 
o"2/o"i and a-s/a^, is expected to be significantly smaller than for hydrogen. 
Ray tracing simulations [8] show that the deviation is in fact negligible. From 
these consideration it follows that the only transmission ratio that needs to 
be measured is r2i = 02! ox for hydrogen. 

From the BRP calibration the transition efficiencies are extracted as a func- 
tion of the parameter r^x- The transmission ratio can thus be evaluated by 
minimizing the of the calibration results. As the transmission probability 
of the sextupole system depends on the velocity distribution of the atomic 
sample, the storage cell temperature T^eii has to stay constant over the cal- 
ibration measurement. This procedure applied to a calibration measurement 
taken at Tceii =95 K has produced the following result [11]: 



02 



= 1.029 ±0.0015. (24) 



The small uncertainty of the measured ratio can be neglected in the calculation 
of the systematic uncertainty of the polarization measurement. The potential 
cell temperature dependence of the transition efficiencies can be neglected as 
the velocity distribution of the sample is fixed by the BRP sextupole system. 
Therefore, a polarization measurement performed collecting at least one signal 
more than the minimum required number (as is always the case) can be used 
to fit the ratio T2\ for various cell temperatures. 

A more detailed description of the BRP calibration and performance can be 
found in Ref. [11]. 



4 Processes occuring inside the Storage Cell 



Inside the storage cell, the atoms of the target may experience several inter- 
actions that have an infiuence on the average target polarization P^- Their 
complete understanding is therefore fundamental for the reduction of the sys- 
tematic uncertainty of the target polarization. In this section an overview of 
these processes will be given, while a quantitative estimate of their contribu- 
tions for the various running periods of the HERMES target can be found in 
Sec. 6. The processes are divided into two categories — recombination pro- 
cesses and spin relaxation processes. 

Recombination and relaxation by wall collisions are surface effects whereas 
spin exchange collisions and resonant interaction of beam and target gas hap- 
pen in the gas storage cell volume. The operating conditions of the target are 
chosen such that these unwanted effects are suppressed as much as possible. It 
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should be noted that spin relaxation can be studied and monitored with good 
statistics by means of the BRP, whereas the residual polarization after recom- 
bination can be studied only using a scattering process like deeply inelastic 
scattering with a much lower rate. 



4-1 Recombination Processes 

The atoms that are injected into the HERMES storage cell may recombine 
to molecules, possibly keeping part of their nuclear polarization. Recombina- 
tion [12] can in principle take place either on the storage cell walls (surface 
recombination) or in the gas phase (volume recombination). At the low gas 
densities of the HERMES target, volume recombination is negligible. On the 
cell walls two types of recombination processes with surface catalysis can be 
distinguished, which are expected to result in different residual nuclear polar- 
ization in the molecules: 

• the Langmuir-Hinselwood (L-H) process, where two adsorbed thermalized 
atoms react with each other; 

• the Eley-Rideal (E-R) process, where an atom from the gas phase reacts 
directly with an adsorbed atom [25]. 

The residual polarization of the surface atoms recombining through the E-R 
mechanism in the cell was measured with the HERMES spectrometer, making 
use of a well known DIS asymmetry. The result was presented in Ref. [14] and 
will be discussed in section 5.2. 

Since the storage cell temperature is usually kept below 140 K, a small amount 
of water originating from the ABS dissociator freezes on the surface and grows 
into an ice layer on top of the Drifilm surface. The presence of water is im- 
portant since water is known to inhibit hydrogen or deuterium recombination. 
Fresh storage cell coatings are hydrophobic, but secondary particles from the 
HERA beam quickly change the chemical structure of the Drifilm surface, 
causing the surface to become less hydrophobic. This was predicted to occur 
in the longitudinal holding field by simulations of acceleration by the beam 
bunch electric fields of positive ions produced by the beam in the target gas. 
Even when not at cyclotron resonance, the ions stochastically gain energy until 
they hit the cell wall with energies of order 100 eV. Surface doses in the range 
lOGy/s can occur, depending on the magnetic holding field and the beam 
current. Furthermore, the downstream portion of the small cell used since De- 
cember 1999 is exposed to synchrotron radiation that is single-scattered from 
collimator CI. Fortunately, it has been observed that during normal operation, 
the quality of the cell surface gradually recovers and the amount of recombi- 
nation in the cell decreases, thanks to the steady accumulation of ice on the 
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surface. 



A detailed description of the recombination studies performed on the HER- 
MES storage cell can be found in Ref. [12]. 

4-2 Spin relaxation 

Three different atomic spin relaxation mechanisms can in principle be observed 
to occur in the HERMES target, as they cause a change in the hyperfine popu- 
lation inside the storage cell. One relaxation process takes place on the wall of 
the storage cell. Another is due to two-body spin-exchange collisions of atoms 
in the gas phase. In the third mechanism the high frequency fields generated 
by the passage of the very short HERA beam bunches can cause a depolar- 
ization in the target atoms under certain conditions [13]. Spin relaxation of 
atomic hydrogen by wall or spin exchange collisions with various types of wall 
coatings has been under study for many years in the context of a hydrogen 
maser for low holding fields [6]. For the HERMES target further studies of 
the magnetic field dependence of spin relaxation as well as studies of the tem- 
perature dependence and density dependence of the transition spectra have 
been carried out [16,30,31,32,33]. The results of these studies are summarized 
in the following sections. 

4-2.1 Wall relaxation 

The atoms injected by the atomic beam source experience several hundred wall 
collisions during their diffusion process through the storage cell. An attractive 
Van-der-Waals potential between the storage cell wall and a diffusing atom 
causes the atom to be physisorbed, i.e. to stay for a certain time close to the 
surface until it is desorbed again. The mean residence time on the wall is 
described by the Arrhenius equation: 



where Tq is the high temperature limit of t^, Ef^ and T are the adsorption energy 
and the temperature respectively, and k is the Boltzmann constant. During 
the time an adsorbed atom stays on the wall (~ 10~^° s for the HERMES 
cell), the magnetic moment of its unpaired electron interacts with magnetic 
moments on the surface by means of a dipolar coupling and/or with other 
unpaired electrons through exchange interactions [18,19,20,21]. Depending on 
the surface type, temperature and magnetic field, one can define a probability 
that the atom leaves the wall in a certain hyperfine state. This is described by 
means of transition probabilities Wjk, each of which represents the probability 




(25) 
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to find an atom in state | k ) that was in state | j ) prior to the wall collision. 
The theoretical description and the calculation of the transition probability 
matrix for the HERMES target can be found in [16]. 



4-2.2 Spin Exchange Depolarization 

Atoms colliding with each other in the gas phase can undergo spin exchange, 
where hyperfine transitions can occur via interaction between the electron 
spins [34]. Due to the conservation of angular momentum, the sum of the 
magnetic quantum numbers mp before and after the spin exchange collision 
has to be the same: 



One can see that if an initial population of the gas consists only of pure state 
1 1 ) or I 3 ) in the case of hydrogen and 1 1 ) or | 4 ) in the case of deuterium, no 
relaxation by means of spin exchange is possible. This is actually not the case 
for the HERMES cell, where to double the density, more than one atomic state 
is injected during normal operation. The effect of spin exchange collisions on 
the hyperfine populations can be described by a tensor Mj^,, which gives the 
probability to find an atom in state | i ) after a collision of two atoms in the 
state I j ) and | k ) prior to the collision: The number of spin exchange collisions 
per unit time and volume is 



where n is the gas density, age is the spin exchange cross section and Vr is 
the relative velocity of the atoms. The density dependence of spin-exchange 
collision relaxation, combined with the knowledge of age [35,36], was used 
to determine the density inside the HERMES cell [16]. The results will be 
reported in Sec. 6. 



4-2.3 Evolution equation for the description of the spin relaxation process 

The various spin relaxations due to wall collisions and spin exchange collisions 
cannot be treated separately. After an interaction with the surface during a 
wall collision, an atom in the gaseous phase can then interact through spin ex- 
change collision and, at the next wall collision, interact again with the surface. 
Since the atoms are constantly exposed to both interactions, an appropriate 
description is given by an evolution equation for the hypefine fractional pop- 



A ^ = . 



(26) 



Nse= n^CTse {Vr ) 



(27) 
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Fig. 11. Measured fractional hyperfine populations vs. magnetic field for hydrogen 
for the various injection modes. The measured values are given by symbols, while 
the lines show a fit using the solutions of eq. (28) 



ulations Ui of the hyperfine state | i ) [16,30] 
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The first two terms of eq. (28) characterize the incoming and outgoing fluxes, 
the third term represents the wall relaxation and the last term the spin ex- 
change relaxation. The effect of recombination is not included, and it is as- 
sumed that the total flux from the ABS is constant with time. In eq. (28), 
three time constants appear: (the diffusion time) is the mean occupation 
time of atoms within the storage cell, r/ is the mean time of flight between two 
wall collisions and Tse is the mean time of flight between two spin exchange 
collisions. For times that are large compared to r^, the relative hyperfine pop- 
ulations converge toward a dynamic equilibrium {steady state). As an example 
of the application of eq. (28), the solution of the system of linear equations in 
the steady state (with ^ = 0) was solved directly in the case of hydrogen or 
with the help of a numerical iteration algorithm in the case of deuterium [16]. 
Fig. 11 represents the measured hyperfine fractional populations as a function 
of the magnetic field for hydrogen and a fit to the data using the solutions of 
eq. (28). The fitted parameters are the three time constants and optionally the 
injected atomic state fractional populations. The R and M transition prob- 
ability matrices are calculated on the basis of theoretical models [16]. These 
data show that a holding field greater than 300 mT is needed to have sizable 
polarization for the case of two injected states. 
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In these references, this equation is designated as the master equation 
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4-2.4 Bunch field induced depolarization 



= -^^-'--^^-^ . (30) 



The beam of the HERA electron storage ring is bunched with time period 
Tbunch = 96.1ns between two adjacent bunches. Typically, most of the 220 RF 
buckets are filled, so that the bunch frequency is given by: 

r^bunch = = 10.42 MHz . (29) 

Tbunch 

A large number of harmonics contribute to the induced magnetic high fre- 
quency field close to the HERA beam, because the bunches are very short 
(cr = 30 ps). The transition frequency between two hyperfine states \i) and 
I j ) is given by: 

E,iB)-_E,iB) 
h 

If matches one of these harmonics at a certain holding field B, resonant 
depolarization occurs. The location of the beam-induced resonances is shown 
in Fig. 12. Bunch field induced depolarization at the HERMES target has been 
studied in the cases of hydrogen [13,33] and deuterium [16] for the longitudi- 
nal orientation of the target holding field (before 2001). During this time the 
observed vr (Am = ±1) transitions were sufficiently well separated so that at 
the working field (335 mT) no beam induced depolarization could be found. 
In the case of the transverse field that has been in use since October 2001, 
both a (Am = 0) and tt (Am = ±1) transitions are allowed because the os- 
cillating field surrounding the HERA beam has components both parallel and 
perpendicular to the static holding field. The spacing between two adjacent 
resonances is narrower than in the longitudinal case (only 0.37 mT difference 
in the magnetic holding field aX B = 0.3 T), so a better homogeneity in the 
holding field is needed (see Fig. 12). 

Thus a transverse field value for which beam resonant depolarization can be 
completely suppressed does not exist. A pair of additional correction coils 
mounted closely around the storage cell to improve the field homogeneity was 
installed for the 2003 running to further reduce resonant depolarization. The 
results will be presented in a forthcoming paper. 



5 Target polarization 



The average polarization Pt of the nucleons in the target is a fundamental 
parameter of the HERMES experiment. The problem is to relate the measured 
atomic polarizations in the sample beam with the density-weighted average 
polarization of the nuclei in the cell. The determination of the target polar- 
ization involves several terms that mainly account for two effects: 
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Fig. 12. Possible beam- induced nuclear depolarizing resonances in the HERMES 
target. The frequency difference between pairs of hyperfine states whose transitions 
lead to nuclear depolarization are plotted as a function of the holding field. The 
frequency values are normalized to vhera = 10.42 MHz. The symbols representing 
the resonance conditions are clearly distinguishable for the vr transitions, while they 
overlap with each other in the case of a and vr | 3 ) ^ | 6 ) transitions, which are sep- 
arated by a difference in the intensity of the holding field AS-p of only 0.37 niT. The 
dashed lines represent the working fields of the transversely (left) and longitudinally 
(right) polarized targets. 



• Presence of molecules in the cell, both unpolarized from the residual gas of 
the scattering chamber and/or directly injected from the ABS, and polarized 
molelcules coming from recombination of polarized atoms on the cell surface. 

• Relaxation and recombination effects in the storage cell leading to an inho- 
mogeneous distribution of the atomic fraction and of the atomic polariza- 
tion. 



31 



The average target polarization Pt as seen by the electron beam is given by 
the following equation: 

= aoarPa + O!o{l- ar)pPa, (31) 

where represents the initial fraction of nucleons in atoms in the gas injected 
by the ABS (or atomic fraction in the absence of recombination), Or denotes 
the fraction of nucleons in atoms that have survived recombination, and (3 = 
Pm/Pa represents the ratio of the nuclear polarization of the molecules Pm 
produced by recombination relative to the nuclear polarization of the atoms 
Pa- The knowledge of /5 for the experimental conditions at the HERMES 
experiment will be discussed in Sec. 5.2. 

The values for oq, ar and Pa are calculated using the measurements of the 
TGA and BRP combined with various calibrations. As both detectors measure 
the properties of a sample of the target gas from the center of the target cell, 
it is necessary to relate the results to a corresponding value averaged along 
the storage cell. This is achieved by evaluating the sampling corrections [15]: 



a. 



r 



(32) 



p 

''p = -BMp- (33) 

a 

The sampling corrections and their systematic uncertainties are functions of 
the measured values of a^^"^ and Pa^^, respectively. The functions depend 
on the geometry of the storage cell, its surface properties and the detectors' 
acceptances. They are calculated with the help of Monte Carlo simulations 
of the stochastic motion of particles in the storage cell, where the histories 
of these particles are stored and analyzed. Another possibility is to calculate 
sampling corrections analytically by describing the diffusion process of the 
target gas inside the storage cell with the use of a one- dimensional diffusion 
equation. The issue was addressed in Ref. [15] and will be discussed here in 
Sec. 5.4. 

The various quantities entering the calculation of the average polarization are 
summarized in Tab. 6. 



5. 1 Determination of the quantities a^. and 



In order to calculate the target polarization following eq. (31), it is necessary 
to separate the atomic fraction that would be seen by the HERA electron 
beam if no cell wall recombination were to occur and the fraction of the 
initially injected atoms that survive wall recombination, as seen by the HERA 
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d.Giisity-s<vGrag6(l nuclear pola.riza,tion in the cell 




atomic fraction in the absence of recombination within the cell 


TGA 


fraction of atoms surviving recombination in the cell 


Cq, 


sampling correction for the atomic fraction 




measured atomic fraction in the sample beam 


— Cpr-a 


Til 1 1 O O T" 1 "\ r^loT'l'VOTli~\Yl r\T OT r^T^T c iir\ Trio /^oll 

liucieai puidi iZ;ciLi(jii oi dLuiiio in une ceii 


CP 


sampling correction for the atomic polarizations 


pBRP 


measured nuclear polarization of atoms in the sample beam 


/3 = Pm/Pa 


relative nuclear polarization of recombined atoms 



Table 6 

Meaning of the quantities contained in eq. (31). 

electron beam. The measured quantity a^^^ can also be factored in two terms 
and can be written as (see eq. (9)): 



a 



TGA 



a. 



TGA^TGA 







a,. 



(34) 



Since dissociative ionization can be neglected [10], the only source of atomic 
flux 0a into the gas analyzer is that originating from the atomic beam source. 
The molecular flux may come from three different sources: 

• Undissociated molecules ballistically injected by the ABS. Since the effi- 
ciency of the dissociator is around 80 %, a significant number of molecules 
pass through the nozzle towards the polarizing stage, affected by neither the 
sextupole magnets nor the transition units. Although the majority of these 
molecules hit the vacuum chamber wall and are pumped away, a small num- 
ber enter the storage cell via the injection tube. These unpolarized molecules 
have a triangular density distribution along the cell axis, with the maximum 
density at the center of the cell. 

• Rest gas in the target chamber. These molecules originate mainly from the 
recombination of atoms that have left the storage cell and have thermal- 
ized on the vacuum chamber walls. This molecular component is clearly 
unpolarized and results in a uniform density component inside the storage 
cell. 

• Atoms that recombine into molecules while in the cell. As explained in 
Sec. 2.3, the cell material and temperature are optimized in order to limit 
the depolarization and recombination effects in the target atomic gas. Never- 
theless, atoms can stick on the cell surface and recombine, and the resulting 
molecules may carry a residual nuclear polarization. 
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The normalized total molecular flow rate 0m of atoms into the TGA ionizer 
can be decomposed in the following way: 

0m = 4>bal + 4>rg + 4>r , (35) 

where (pbai, 4'rg and (pr represent the ballistic, the rest gas and the recombina- 
tion components, respectively. To separate between polarized and unpolarized 
molecules, we rewrite the two factors in eq. (34) using the three contributions 
appearing in eq. (35): 



«r = — — — ' (36) 

0a + 0r 



° 0a + 0r + 4>bal + (f^rg 



The independent measurement of each of the three components of the total 
molecular flow rate is accomplished by means of routine target calibrations, 
which are usually performed between two HERA electron fills when the HER- 
MES spectrometer is not taking data. The rest gas contribution 0,.^ to the 
molecular flux is linearly proportional to the pressure in the target chamber. 
The constant of proportionality is determined by injecting into the target 
chamber increasing quantities of hydrogen (deuterium) molecules with a con- 
trolled flow system and by measuring the corresponding TGA count rates. The 
ballistic flux contribution (pbai can be measured by varying the atomic flux of 
the injected beam while keeping the molecular flux constant. The atomic flux 
is changed by using the hyperfine transitions to select either 1 or 2 (or 3 in the 
case of deuterium) hyperfine states. Since the recombination probability in- 
side the cell was found to be independent of the injected beam intensity [12], 
the quantity (pbai can be extracted. Eq. (35) is then solved for the remain- 
ing molecular flux 0^,. After the molecular contributions 0^^, (pbai and 0^ are 
extracted, one can determine a^^^ using eq. (36) and thus a,, using eq. (32) 
with the sampling correction Cq, determined as described in Sec. 5.4. A detailed 
description of how these calibrations are performed was given in Ref. [10]. 

The factor is defined in terms of the nucleon thicknesses ta,r,bai,rg in the cell 
seen by the HERA beam: 



ta tr thai trg 

It can be determined without reference to additional sampling corrections. All 
the 0j terms in eq. (37) represent nucleon fluxes. Let us define a proportionality 
constant c representing the TGA response to atoms through the relationship 
0a = eta- The fluxes (pr, 4>bai and 0^^ are suppressed by a factor of \/2 due to 
the smaller molecular velocity, leading to 0^ = ctrj^f^ and 0fea; = ctbai/V^- 
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Furthermore, in the hmiting case of small cell wall recombination, the densities 
of atoms and molecules producing ta, thai and tr are triangular. In contrast, 
the rest gas producing trg has a uniform density along the cell, resulting in 
= ctrg /{2\/2). If the difference between the conductances of the atoms and 
molecules is neglected so that they have the same distributions along the cell, 
the above considerations can be combined to give 

(pa + V2<Pr 

ao = 1=^ . (39) 

(f)a + V2(0r + (j)bal + 20^^) 

However, this difference in conductances can be taken into account [16,24], 
leading to different distributions along the cell and a more complicated ex- 
pression: 

di(f>a + <Pr) 

Oo = ^ 7= 5 (,4Uj 

d{(t)a + 0r) + y'2{(pbal + 20r.g) 



where d = ^2/(1 + {^2-1] 



Or 



5.2 Relative polarization (3 of nuclei in recomhined atoms 



As mentioned in the previous sections, some of the polarized hydrogen or 
deuterium atoms may recombine in the storage cell. A direct measurement of 
the remnant polarization contained in the molecules is not possible at HER- 
MES, as the Breit-Rabi polarimeter can measure the atomic polarizations 
only. Atomic recombination is a surface-mediated temperature-dependent pro- 
cess [12,16]. Hence measurements of the nuclear polarization of recombined 
atoms using storage cells of various materials and/or with various coatings 
[37,38] are not directly applicable to the HERMES case. Without any infor- 
mation about the polarization (3 of nucleons in molecules (or molecular polar- 
ization in short) one has to allow in eq. (31) the conservative range < /3 < 1. 
A set of data taken in 1997 at a higher cell temperature (260 K instead of 
the nominal 100 K) was used to measure j3 for hydrogen at that temperature, 
resulting in /^^eoi^ = 0.68±0.09,tat ±0.06,3/,t [14]. This result can be used as an 
upper limit for (3 at 100 K {l3l^g,f) under the following assumptions. The main 
mechanism responsible for recombination at 100 K and 260 K is the same (i.e. 
the FrK mechanism) [12]. At 100 K, the recombination probability is smaller 
than at 260 K, because the atoms impinging from the volume have less kinetic 
energy to overcome the activation barrier. The residence time of chemically 
adsorbed atoms on the surface at 100 K is greater than at 260 K, and their 
possible polarization cannot exceed the value measured at 260 K. This results 
in = 0.83. 

The lower limit for /3^^^^ iPiow^) can be derived from a simple argument. By 
assuming that the nucleon spins are not affected by the recombination process. 
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the nuclear polarization of the molecule at the time of its formation (P^) can 
be evaluated by taking the average value of the polarization of an atom coming 
from the volume (Pa) and one resident on the surface (Ps)- The theoretical 
lower limit for f3^^^^ can then be inferred by supposing that atoms sticking 
on the cell surface are totally depolarized (p^ioo^''°"' = Q) and applying the 
equation: 



plOOK,low 
alOOK _ m 
Plow 



plOOK 
plOOK _|_ plOOK,low 2 



plOOK 



(41) 



The result is Plo^^ = 0.5. A factor of 0.9 takes into account the depolarization 
of molecules while colliding with the walls [14,38], ending in (3}^^^ = 0.45. 
Hence, the range of (3^^^^ can be limited to 

0.45 < < 0.83 . (42) 

The corresponding value and uncertainty of j3^^^^ are therefore given by: 

/^looi^ = 0.64 ±0.19. (43) 



5.3 Injected atomic polarizations 



An important parameter for the evaluation of the target polarization is the 
injected atomic polarization P™-' (i.e. the nuclear polarization of atoms in the 
absence of depolarization), as it enables the separation of the various spin- 
relaxation effects. The injected polarization for each ABS operating mode 
(Tab. 2) can be calculated if the transmission probabilities of the sextupole 
system and the efficiencies of the adiabatic transition units are known. The 
sextupole transmission probabilities for hydrogen and deuterium have been 
calculated using a Monte Carlo simulation and the results are given in Ref. [8]. 
The injected polarization is fairly insensitive to these values. The efficiencies 
of the transitions of the ABS can be extracted by measuring the magnetic 
field dependence of the spin relaxation processes and by fitting the data with 
a theoretical model described in Sec. 4.2.3 and given by eq. (28). For normal 
operation, the important ABS injection modes are those having a large positive 
or negative nuclear polarization, Pz+ and respectively, and, in the case of 
deuterium, the two additional positive and negative tensor polarization {Pzz+ 
and Pzz-) modes (see Tab. 2). 

For hydrogen, only two transitions after the sextupole system are used: the 
SFT 2-4 for P^_|_ and WFT 1-3 for P^_. Therefore, the injected polarizations are 
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directly given by the transition efficiencies whose values have been extracted 
in Ref. [16] with the technique mentioned above, and cross-checked using a 
more direct method based on the symmetry of the spin relaxation processes 
for the two injection modes. In the case of deuterium, several transition units 
are used to provide the four polarizations, some of which are located between 
the two sextupole subsystems. Hence, the injected polarizations have to be 
calculated taking into account the transition efficiencies as well as the trans- 
mission probabilities for atoms whose hyperfine state gets changed between 
the two subsystems. 

The calculated polarization for the hydrogen and deuterium beams injected 
into the target are summarized in Tab. 7. The degree of nuclear polarization 
injected by the ABS when running with deuterium is substantially lower than 
for hydrogen. The reason is that, due to the higher number of hyperfine states, 
the deuterium nuclear polarization can be achieved only by using the transition 
units located between the two magnetic subsystems, resulting in an incomplete 
separation of states. 



5.4 Sampling Corrections 

The sampling correction functions Ca{c(^'^^) and cp{P^^^) depend strongly 
on the distribution of recombination and wall relaxation probabilities per wall 
collision along the storage cell. For a freshly produced storage cell that has not 
yet been exposed to the HERA beam, one can assume that the recombination 
and wall relaxation probabilities are constant over the entire surface of the 
storage cell, or at least that microscopic inhomogeneities due to various types 
of surface sites are equally distributed. Possible mechanisms by which the 
HERA beam could affect the surface properties of the storage cell include: 



Mode 


Hydrogen 


Deuterium 


pinj 


97.26 ±0.05 


87.78 ± 0.07 


^z — 


-97.38 ±0.07 


-89.50 ±0.10 






96.22 ±0.38 


pinj 

zz~ 




-177.30 ±0.61 



Table 7 

Estimated polarizations for the hydrogen and deuterium atomic beams injected by 
the ABS into the target. All values are multiplied by 100. The quoted uncertainties 
take only the dominant uncertainties of the ABS efficiencies into account. The values 
refer to the polarization of injected gas with the 0.3 T holding field. Note that, 
due to the standard way of defining it, the maximum achievable negative tensor 
polarization is -200. 
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• desorption of the surface coating by synchrotron radiation; 

• bombardment of the storage cell walls with charged particles leading to 
radiation damage; 

• chemical reactions of ionized hydrogen or deuterium radicals with the cell 
surface; 

• redistribution of the ice layer on the storage cell by RF heating due to the 
wake fields; 

• contamination of the storage cell surface by sputtered material from the 
collimator upstream of the target. 

Under normal operating conditions no influence of the HERA beam on the 
target performance has been observed. Nevertheless certain distinct incidents 
affecting the cell have been observed, when the beam tune was not optimal 
or the HERA beam was accidentally lost near the HERMES target region. 
Since the effect of the HERA beam on the surface properties is poorly known, 
every realistic distribution of recombination and wall relaxation probabilities 
has to be considered when evaluating the sampling corrections. It is therefore 
convenient to distinguish certain scenarios of special interest. 

• Homogeneous cell (HC): the natural assumption for a freshly produced 
cell not yet exposed to the HERA beam. The calculated sampling corrections 
set an upper limit for a,, and Pa for any given values of aj'^^ and P^^^ , 
respectively. 

• Homogeneous beam tube (HBT): the main cell tube is homogeneously 
affected by the beam, while the side tubes are unaffected. 

• Inhomogeneous beam tube (IBT): The main cell tube is inhomoge- 
neously affected by the beam, while the side tubes are unchanged. This 
scenario does not lead to a single sampling correction, but it can be used 
to calculate a lower limit for or Pa for any given value of aj"^^ or P^^^ , 
respectively. 

A model for calculating the sampling corrections was developed in Ref. [16]. 
It is based on a parameter 7r(7d), the recombination (spin- flip) probability 
during a wall collision. If these probabilities were constant and equal for the 
beam and sample tubes, the sampling corrections would depend only on their 
geometry. This is actually the case for a freshly installed cell, whose Drifilm 
coating is known to be intact and uniform. In all scenarios it is assumed that 
the relation 7^^-^ < 7^-^" holds at any time for the recombination probabilities 
in the sample tube, •y^'^, and beam tube, 7^"^. The validity of this hypothesis is 
proven by the data [16,24]. The numerical calculations of the sampling correc- 
tions have been performed using a molecular flow simulation, complemented 
by analytic calculations based on the one- dimensional diffusion equation [16]. 

The sampling correction cp must be subdivided into two terms, one for the 
wall collisions and a second term for spin exchange relaxation, as the two 
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processes must be treated in different ways. This can be done by defining the 
depolarization factors (wall collisions) and D^e (spin exchange relaxation) 
so that for small spin relaxation, one has: 

CP = " ~ i^^rf Use 

^ pBRP r)BRP DBRP p^-^i ' ^ ' 

As the spin exchange relaxation occurs in the gas phase, and for the typical 
values of > 0.9 the atomic density distribution is little affected by any 
change of the cell surface properties, the relation between Dse and DfJ^^ is 
similar for all of the listed scenarios. Moreover, since calculations have shown 
that — Df^^ [16], one can conclude that the sampling correction cp is 
determined by the wall depolarization only: 

Pa Dwd f Ar\ 

= pBRP - nBRP ■ (45) 

-'a ^wd 

As an example, in Fig. 13 the effect of the calculated sampling corrections for 
the large-aperture cell used in the 1996-97 running with hydrogen is shown. 
During normal operation, the values of aj'^"^ and D^f^ are found to be almost 
always well above 0.9, so that local linear approximations may be fitted to the 
calculated sampling corrections: 



ar = a+{l-a) aj^^ (46) 

= a^«^+a(l-a^^^), (47) 

D^, = 6+(l-6)Z^fr (48) 

= D^r+bil-D^r)- (49) 

The values of the parameters a for recombination and b for depolarization are 
fitted to the average between the two extreme scenarios shown in Fig. 13, and 
their systematic uncertainties Ad and Ab are taken to be the differences found 
when fitting the extremes. 



5. 5 Evaluation of the systematic uncertainty in Pt 



After having determined all the terms of eq. (31), Pt can be expressed as 
a function of the measured quantities. It is convenient to write the target 
polarization in the following way: 

Pt = ao {ar + (1 - a,-) P)Pa = oieff ■ Pa ■ (50) 
where the term a^ff is called effective atomic fraction. 
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Fig. 13. Simulated correlation between the values measured by the TGA or BRP 
detectors and their corresponding density-weighted average values in the storage 
cell. The three solid curves represent the three scenarios for the state of the cell 
surface: homogeneous cell (HC), homogeneous beam tube (HBT), and inhomoge- 
neous beam tube (IBT). The left figure shows the range for recombination and the 
right picture shows the range for wall depolarization. Acceptable uncertainties of 
Or and Dwjj (differences among the three horizontal dotted lines) can be achieved 
only with aj'^^ and D^j^ close to unity (vertical dotted lines). The corresponding 
uncertainty grows rapidly with increasing recombination or depolarization. 

5.5.1 Uncertainty of the effective atomic fraction a^ff 

Combining eqs. (50), (47), (40), (36) and (35), tte// can be expressed as 

_ 4>a + (g (1 - /3) + P){(pm. - 4>bal - <Prg) 
"'^^ ~ (l)a + K + {^- l)0bai + (2^2 - 1)0., ■ 

The systematic uncertainty of ae// is determined by the measured flow rates 
0a; 0m; 'Phai ^ud 0rg and their uncertainties as well as the values and uncertain- 
ties of (3, a and the TGA calibration constant k. With the matrix T deflned 
by 

--(^•••••%^)' 

and the diagonal covariance matrix C of all calibration constants, the system- 
atic uncertainty in Oe// is given by 

Asysaeff = VTCT^ . (53) 

The matrix T contains 10 independent derivatives, i.e. for k, a, (3 and the 
7 additional terms related to the ballistic flux and rest gas calibrations as 
described in Ref. [10]. 
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5.5.2 Uncertainty of the atomic polarizations Pa 

Combining eqs. (44), (45) and (49), the average nuclear polarization Pa of 
atoms in the HERA beam can be written as 

Pa = Pa^'''' + b {Dse P^ ' Pa'''') (54) 

For the calculation of the systematic uncertainty Asys Pa, the important un- 
certainty sources to be included are Ab and AgygP^^^^^ . As a result, one 
finds: 

^sysPa = {{l-m^sysPa^n'+il-D^^dfDliPrrAb'y^". (55) 

During the analysis one typically averages over a large dataset, so that the 
statistical uncertainty of the polarization measurement Agta Pa is negligible. 



6 Performance 

In the present section an overview is given about the target performance in 
the running period 1996-2003, by analyzing distinct periods for each target: 

• 1997 for the longitudinal hydrogen target; 

• 2000 for the longitudinal deuterium target; 

• 2002 for the transverse hydrogen target. 

6.1 Longitudinal hydrogen running in 1997 

The upper plot in Fig. 14 shows the atomic fraction a^'^^ measured by the 
target gas analyzer during the 1997 running period. The nuclear polarization 
of atoms measured by the BRP is plotted in the lower half of Fig. 14. Aside 
from startup problems, the target performance was relatively smooth. Unfor- 
tunately, a severe HERA beam loss close to the HERMES region occurred at 
a certain point (left line), resulting in a change of the cell surface properties. 
After this bad event, it was decided to collect a set of data at higher tem- 
perature {Tceii = 260 K instead of the nominal value Tceii = 100 K) in order 
to measure the (3 parameter (see Ref. [14]). At the end of this period, the 
storage cell was exchanged (middle line). Also the second cell experienced a 
beam loss (right line) which produced a small hole in the cell wall. For these 
reasons, the analysis of the 1997 dataset was divided into four periods, during 
each of which the cell surface condition was assumed to be constant. For each 

The uncertainties on APge and on Pl^-' are negligible. 
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Fig. 14. The TGA (above) and BRP (below) measurements for the entire 1997 data 
taking period. The vertical dashed lines indicate HERA beam loss events which 
affected the cell surface properties, and the solid line indicates the replacement of 
the target cell. 



period, specific studies have been performed to limit the uncertainties on 
and Pa [24]. 

For the entire running period with hydrogen, we quote a luminosity-weighted 
average value for the target polarization of = —Pz- = 0.852 ±0.033 and a 
density of 7.6 x 10^^ nucl./cm^ . The latter number was derived from a method 
exploiting the density dependence of spin exchange collisions [16]. 
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6.2 Longitudinal deuterium running in 2000 



Over the eight month period of deuterium running in 2000, the HERMES tar- 
get operated under very stable conditions. The atomic beam source, equipped 
with a microwave dissociator, experienced only one major failure, which hardly 
affected data accumulation. Neither the scattering chamber nor the TGA/BRP 
vacuum system had to be opened, leading to constant amounts of residual gas 
in the storage cell as well as stable efficiencies of the detectors monitoring the 
target. There was no need to replace the storage cell, as its performance in 
terms of atomic recombination and nuclear depolarization was always excel- 
lent. The reason for such good behavior can be found both in a particularly 
fortunate Drifilm coating and in the limited damage caused by a very well 
tuned electron beam, i.e. no unwanted beam losses in the vicinity of the tar- 
get. 

The outstanding performance of the HERMES target in this period is well rep- 
resented in Fig. 15 (note the change in scale compared to Fig. 14). The stable 
behavior of the cell surface in terms of atomic recombination is demonstrated 
in the upper plot of Fig. 15, where the atomic fraction measured by the target 
gas analyzer is plotted. The lower plot shows the nuclear polarization. In the 
period between January and June 2000, the target was operated with only 
positive and negative states of nuclear vector polarization, implying positive 
tensor polarization. From July on, a new injection mode combining vector 
and tensor polarization Pz+, Pz-, Pzz+ and Pzz- was established in order to 
perform the first measurement of the tensor-polarized structure function bf. 

The average values for the four polarization states used in 2000 are the fol- 
lowing: 



P-. 



z+ ■ 



p. 



ZZ+ ■ 



-0.851 ±0.029 
-0.840 ± 0.026 
-0.891 ±0.027 
-1.656 ±0.049, 



(56) 
(57) 
(58) 
(59) 



Note that only the systematic uncertainty is quoted as the statistical un- 
certainty is negligible in comparison. The differences between the systematic 
uncertainties originate from the injection mode dependent systematic uncer- 
tainty of the atomic polarizations. The measured target density for this pe- 
riod, measured by means of the spin-exchange collision relaxation method, is 
2.1 X 10^^ nucL/cm^. 
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Fig. 15. Atomic fraction, a^^"^ upper plot, and vector/tensor polarization, p^^P 
lower plots, measured by TGA and BRP respectively during the year 2000 running 
period. Each symbol represents data averaged over a 72-hour bin. The absolute 
bin averaged values are always above 0.9. The tensor polarization (lower plot) was 
employed from July on. 

6.2.1 Spin relaxation 

The spin-flip probability per wall collision 7^ is proportional to r^, the mean 

residence time: Tg = tqc^^'^' (see eq. (25)). Therefore 7^ depends on the wall 
temperature T. In other words, if P^^^ does not change with T, wall depolar- 
ization effects are vanishingly small within the experimental uncertainties. In 
Fig. 16, a scan of the cell temperature down to very low temperatures is shown. 
The result indicates that no wall depolarization occurs. As a consequence, for 
the deuterium target the sampling correction for nuclear depolarization equals 
unity (cp = 1), so that: 

Pa = P^''''. (60) 
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Fig. 16. Nuclear vector polarization p^^^ as function of the cell temperature T 
measured in August 2000. No dependence on the temperature is observed. The 
operating temperature of the target during data collection was 60 K. 

Moreover, the calculated depolarization due to spin exchange collisions re- 
sulted in AP^'^ = 0.3 %. As in 2000 spin relaxation can be neglected com- 
pletely, the polarization injected by the ABS equals the measured BRP value: 

pinj ^ pBRP _ ^g^^) 

The polarizations measured by the BRP were found to be consistent with 
those calculated as described in Sec. 5.3, within the uncertainties. 



6.2.2 Recombination 

Although atomic recombination depends strongly on temperature, no variation 
of the measured atomic fraction a^*^"^ over a very wide temperature range 
around the working temperature Tceii = 60 K was ever observed during the 
whole year. Fig. 17 (upper plot) reports a cell temperature scan performed 
in August. (For comparison a temperature scan performed with hydrogen in 
1997 is also shown in the lower plot of the same figure). Additionally the 
measurement of the fraction aj'^^ of atoms surviving recombination yielded 
the following average value: 

a^^^ = 0.997 ±0.0U, (62) 

confirming that within the measured uncertainty no evidence of recombination 
was found. Therefore, for the deuterium operating period in 2000, both the 
deviation from unity of the sampling correction Cq, and the impact of the 
uncertainty in jS are negligible. 



6.3 Transverse hydrogen running in 2002 



During the 2002-2003 operating period, the HERMES transversely polarized 
hydrogen target worked in a very stable way. On the other hand, the poor 
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Fig. 17. Upper plot: measured atomic fraction as a function of the cell temperature 
for Tceii < 100 K for the deuterium running in 2000. No dependence on the tem- 
perature can be seen. The working temperature during normal operation was set 
to Tceii = 60 K. For comparison a temperature scan with hydrogen taken in 1997 is 
shown in the lower plot. 

performance of the HERA beam did not allow a detailed study of the bunch 
induced depolarization. Nevertheless, by scanning the target magnetic field a 
working field value was established where the loss of polarization due to the 
field non-uniformity was limited to 1.5 % at the highest available beam currents 
(25 mA). The average target polarization {Pz+ — Pz-)/^ during the 2002 data 
taking period was 0.783 ±0.041. This value is lower than that reported for the 
longitudinal hydrogen target (Section 6.1), mainly due to the lower applied 
target holding field and to the fact that the higher density allowed by the use 
of the intermediate size cell resulted in an increase of spin relaxation due to 
spin exchange collisions. 



The stability of the transverse target operation during the 2002 running period 
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Fig. 18. The TGA (above) and BRP (below) measurements for the 2002/03 data 
taking period. 

is illustrated in Fig. 18. The target density in this period measured via spin- 
exchange collisions was 1.1 x 10^^ nucl./cm^. 



7 Discussion 



Tab. 8 summarizes averaged values for the target properties during the various 
running periods. A comparison between the longitudinally polarized hydrogen 
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Table 8 

Summary of the HERMES target performance in the years 1997, 2000 and 2002. 
The tabulated values are averaged over the respective running periods. The APj's 
represent the polarization losses arising from the listed process (see text). 



and deuterium target leads to the following conclusions. At the same guide 
field, in the case of deuterium the spin exchange and wall relaxation processes 
are suppressed by a factor [B^jB^Y — 20 compared to hydrogen. The posi- 
tive and negative Pz values for hydrogen coincide very well, while for deuterium 
there is a clear difference. This is due to the larger number of efficiencies and 
transmissions involved. 

From the table it is possible to deduce departures from the optimal perfor- 
mance like the amount of recombination (Aa^ = 1 — a^), and the depolar- 
ization due to the various depolarizing mechansims: spin-exchange collisions 
(APse), wall-depolarization (APwd) and beam-induced fields (APb/)- The 
target density integrated over the cell length, or areal density t, is also shown 
together with a relative factor of merit (FOM) that takes into account both 
the polarization and the density. It is worthwhile to mention again that in 
1997 the target cell had a larger cross-section (29mmx9.8mm) implying a 
higher conductance and a correspondingly lower density (about 40%). The 
table also illustrates once more the optimal performance of the D\\ target in 
2000. For this target no recombination (Aa^ = 0) and no depolarization phe- 
nomena {APw£) = APsE = APbi = 0) have been observed. This excellent 
behavior has to be attributed to the lower critical field of deuterium (11.7 mT) 
compared to hydrogen (50.7 mT), and to an improvement in the knowledge 
and preparation of the cell surface coating. The improved surface coating is 
also responsible for the better performance of the 2002 H± target compared 
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Table 9 

Deeply inelastic scattering events collected with the three polarized targets after 
data quality selection. The numbers are given in millions. 

to the 1997 H\\ target, as far as the amount of recombination is concerned. 

On the other hand, in 2002, due the higher density obtained and the higher 
number of wall collisions, spin- relaxation effects became more important, caus- 
ing a loss in polarization of about 11 % and limiting the gain in the factor of 
merit with respect to the 1997 running due to the increase in density. This 
indicates how an enhancement in the figure of merit of the target might be 
obtained by increasing the holding field proportionally with the density, in 
order to limit the spin relaxation effects. 



8 Selected HERMES results 

In this last Section before the conclusions, some of the measurements which 
have been carried out at HERMES with the polarized internal gaseous target 
during the various running periods are mentioned. The number of deeply in- 
elastic scattering events (after data quality selection) collected with the three 
polarized targets is given in Tab. 9. 

Among the measurements which have been performed using the HERMES 
detector, we quote the following few: the measurement of the Qi structure 
function [39]; the measurement of the gf structure function [40]; the first mea- 
surement of the bf structure function [41] ; measurements of the quark-fiavour 
decomposition of the nucleon spin [42,43]; the first measurements of single-spin 
asymmetries in deep-inelastic scattering off both longitudinal and transversely 
polarized targets [44,45,46] and a first study of the gluon polarization in nu- 
cleons [47]. 



9 Conclusion 

The HERMES polarized hydrogen/deuterium target has been reliably operat- 
ing since 1996. It provides polarization values between 0.78 and 0.85, densities 
1-2 X 10^^ nucleons/cm^, and total systematic uncertainties of only 0.03-0.04. 
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The systematics of both the recombination and depolarization processes in 
the storage cell are fully under control. 
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